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We report a Raman scattering investigation of multiferroic bismuth ferrite (BiFeOa) epitaxial 
(c-axis oriented) thin films from -192 to 1000°C. Phonon anomalies have been observed in three 
temperature regions: in the 7-phase from 930°C to 950°C; at ~ 370°C, Neel temperature (Tn), and 
at ~ -123°C, due to a phase transition of unknown type (magnetic or structural). An attempt has 
been made to understand the origin of the weak phonon-magnon coupling and the dynamics of the 
phase sequence. The disappearance of several Raman modes at ~ 820° C (Tc) is compatible with 
the known structural phase transition and the Pbnm orthoferrite space group assigned by Arnold 
et al. P]. The spectra also revealed a non-cubic /3-phase from 820-930° C and the same non-cubic 
phase extends through the 7-phase between 930-950° C, in agreement with Arnold et al. Qi and an 
evidence of a cubic (5-phase around 1000° C in thin films that is not stable in powder and bulk. Such 
a cubic phase has been theoretically predicted in 01 . Micro- Raman scattering and X-ray diffraction 
showed no structural decomposition in thin films during the thermal cycling from 22-1000° C. 

PACS numbers: 77.55. +f, 78.30.-j, 77.80.Bh, 78.66.Nk 
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I. INTRODUCTION 

Multiferroics are the materials which display a coexis- 
tence of at least two of the switchable states: polariza- 
tion, magnetization or strain in the same phase In 
addition, they may also exhibit a magnetoelectric (ME) 
effect: magnetization induced by an electric field and 
electric polarization by means of magnetic field Q . The 
current interest in multiferroics is largely based on engi- 
neered epitaxial and heterostructured thin films, because 
their physical properties are as good as bulk and per- 
mit technological applications in data storage, magnetic 
recording, spintronics, quantum electromagnets, and sen- 
sors [6|-(8j. Devices made up of multiferroic materials 
can perform more than one task and facilitate device 
miniaturization. A weak ME effect has been observed 
in most multiferroics, generally showing a small change 
in their spontaneous polarization under applied magnetic 
field 1,1 [H. However, the complete switching of fer- 
roelectric domains by applied magnetic fields has rarely 
been observed. Why and under what circumstances a 
large coupling should exist and how to control the cou- 
pling are still open questions. Understanding the physics 
of the different possible interactions between magnetic 
and electric order parameters i.e. giving rise to magne- 
toelectric (ME) coupling would be very useful. 

Magnetism and ferroelectricity are involved with lo- 
cal spin ordering and off-center structural distortions, re- 
spectively . These are quite complementary phenom- 
ena that coexist in certain multiferroic materials. Cur- 
rently, BiFe03(BF0) is one of the most widely studied 
multiferroics because it is one of only two or three single- 
phase multiferroics at room temperature i.e an antifer- 
romagnetic (AFM) incommensurate phase with cycloidal 
modulation (A w 60 nm) below w 370 °C fer- 



roelectric up to ~ 820°C [lj|, and ferroelastic between 
820-930°C {3. Bulk BFO crystalhzes in a rhombohe- 
dral (a = 5.58 A and a = 89.5°) structure at room 
temperature (RT) with space g roup R3c (C3„) and an- 
tiferromagnetism of G-type [ij, [l^ . The structure and 
proper ties of bulk BFO have been studied extensively 
[13. [la - flq and although early values of polarization were 
low (Pr = 6.1 jjiC/cxs?) due to sample quality, Pr = 40- 
100 /iC/cm^ was recently found in bulk by several dif- 
ferent groups [H, [l^l . The epitaxially grown thin films 
of BFO on SrTiOa (STO) substrates show very high 
values of Pr('^ 100 /iC/cm^) Q compared to the best 
known ferroelectrics such as PbZrTiOs (^ 70/iC/cm^) 
and BaSrTiOs (- iQ^J.G|cTc?). This makes BFO a po- 
tential material for novel device applications. 

The motivation for the present study is manifold. The 
first objective is to test the recent space group deter- 
mination of the 7-phase reported by Arnold et al. [l[ 
as being orthorhombic. Their definitive neutron study 
showed that the 7-phase is indeed stable (which in itself 
had been controversial), and that it has the same or- 
thorhombic Pbnm orthoferrite symmetry as does the /3- 
phase. A cubic Pm3m perovskite structure was definitely 
ruled out, although a body- centered orthorhombic space 
group was indistinguishable from the primitive Pbnm. A 
main aim of our Raman study is to test the orthorhom- 
bic crystal class for the 7-phase and see whether we can 
further distinguish between primitive and body-centered 
orthorhombic structures. 



II. EXPERIMENTAL DETAILS 

We investigated 300 nm (001) BFO thin films on STO 
(100) substrates with ~ 25 nm thick SrRuOs (SRO) 
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buffer layer by pulsed laser deposition (PLD). A Jovin 
Yvon T64000 micro-Raman microprobe system with Ar 
ion laser (A = 514.5 nm) in backscattering geometry was 
used for polarized and temperature dependent Raman 
scattering. Sample deposition and experimental details 
are given in [l^ . 

The X-ray diffraction (XRD) pattern (see Fig. EK) of 
the BFO films taken using CuKq, (1.5406A) radiation 
show c-axis (pseudo-cubic <001> direction perpendicu- 
lar to the substrate) orientation with a high degree of 
crystallinity. The c-axis length was found to be 3.95 A, 
which implies epitaxial strain is quite relaxed. This 
agrees with the reported values (c = 3.997 A) [sTj . 

The comparison of the unpolarized (perpendicular to 
the <001> of the substrate) Raman spectrum of BFO 
thin film with STO and SRO/STO spectra (cf. Fig. 2a 
in [l3|) precludes any Raman contribution from the sub- 
strate and bottom electrode; to the contrary, we observed 
a dip, rather than a peak, at the STO strongest peak po- 
sition. As is evident from the intensity comparison, all 
of these peaks are due to the BFO normal modes of vi- 
brations and none of them arose from the substrate. We 
verified our results using target materials, single crys- 
tals, and also by growing (001) BFO films on different 
substrates. 



III. RESULTS AND DISCUSSION 
A. Phonons in the 7-phase 

Fig H] shows the temperature dependent Raman spec- 
tra of (OOl)BFO film (aOOnm thick) on SrTiOa substrates 
with SrRuOs buffer layer (25nm thick). As can be seen 
from Fig 1T]d, four phonon features marked with arrows 
persist into the 7-phase from 930-950°C. Therefore this 
confirms an existence of non-cubic 7-phase. No first- 
order Raman scattering from phonons is allowed for cu- 
bic PmSm (because each ion is at an inversion center, all 
phonons are odd-parity). This is in agreement with the 
observation of an orthorhombic symmetry for the 7-phase 
by Arnold et al. [l| . Although the Raman lines are broad 
and weak at these temperatures, they exhibit no signifi- 
cant frequency shift at the ^-7 transition at 930° C; nor 
are there additional lines above or below 930°C. There- 
fore it is very likely that the structure is Pbnm in both 
/? and 7 phases and that no primitive-to-body-centered 
phase change occurs. Note that Fig lT]is for a thin film. 
Therefore we conclude that the films do not differ from 
single crystals and powders previously studied with re- 
gard to the 7-phase. 

There is renewed interest in this /3-7 transition because 
of the earlier discovery of a metal-insulator transition 
at 931°C where the orthorhombic-cubic transition takes 
place in bulk [isj and that the same transition may oc- 
cur at 47 GPa at RT. Earlier Mossbauer studies estab- 
lished that at this pressure the magnetization also dis- 
appears [211]. Although Gavriliuk et al. [2l| concluded 



that this is a rhombohedral-rhombohedral symmetry- 
preserving Mott transition, that seems quite unlikely, be- 
cause Haumont et al. [l^] have found several phase tran- 
sitions at lower pressure. Thus the symmetry of BiFeOa 
is not rhombohedral on either side of the high-pressure 
metal-insulator transition. Whether the transition is a 
Mott transition or a band transition is unproven. Various 
theoretical models disagree: Vasquez et al. Q get an ah 
initio Mott transition; Clark et al. [2^ got a band tran- 
sition to a semimetal from a screened exchange model. 

Figs, [iji and b show the temperature variation (from 
RT up to 1000°C) of unpolarized Raman spectra of a 
BFO (001) thin film. A closer observation near the phase 
transitions reveals two noticeable changes in the signa- 
ture of the Raman spectra: the disappearance of several 
stronger modes at ~ 820 °C and the complete disappear- 
ance of all the modes above 1000 °C. This temperature 
behavior implies that BFO maintains its room temper- 
ature structure up to ~ 820 °C, indicating the a struc- 
tural (ferroelectric) phase transition, in agreement with 
the earlier investigations on BFO bulk single crystal and 
polycrystalline samples [l^] . Note that thin films of BFO 
show first order phase transitions as in bulk, whereas 
STO and PbTi03 PTO are known to be first order in 
bulk but second order in thin films [2^, [2^ . 

The presence of the four peaks {'^ 213, 272, 820 and 
918 cm-i) above 820°C up to ~ 950°C (Fig. [lb) shows 
that the intermediate beta-phase is not cubic (Pm3m) as 
reported by Haumont et al 2J|. In fact, the phase dia- 
gram of BFO [23, and its more recent revised versions 
|15| . show that BFO possesses a non-cubic /3- phase be- 
tween 820 to 933° C before it goes to the 7-phase, and the 
/3-phase was recently shown to be orthorhombic by using 
high temperature X-ray diffraction and domain struc- 
tures [iBl and neutron diffraction [Ij. The complete dis- 
appearance of peaks at above 950°C - not at 93CP C - 
indicates that the 7 high temperature phase also can- 
not be cubic (PmSm), for which any first-order Raman 
scattering is forbidden. 

Fig. [2^ shows the room-temperature XRD patterns of 
an as-grown film and film after it underwent 1000°C ther- 
mal cycle. As can be seen, the as-grown film is highly 
epitaxial showing only (001) peaks and became polycrys- 
talline after thermal cycling. In principle it is possible 
that the specimen would melt at high temperatures and 
then recrystallize in the specimen holder (bottom of the 
Pt crucible). However, we monitored the sample surface 
continuously with an optical microscope and no thermal 
decomposition was observed up to 1000°C. 

The Raman spectra (Fig. ^p) before and after heat- 
ing show exactly same number of phonon modes, indicat- 
ing either no decomposition up to 1000°C or complete 
recrystallization, contrary to earlier studies [T7| . which 
could be due to the reduced surface/ volume ratio, mini- 
mal surface imperfections and defects, and increased sta- 
bilization from the substrate. Note that the possibili- 
ties of subtle structural changes (small changes on angles 
and/or in-plane lattice parameters) cannot be completely 
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FIG. 1: Temperature dependent Raman spectra of (OOl)BFO film on SRO/STO from 22-400° C (a) and 550-1000° C (b). The graphs 
in (b) were adapted from our earlier work [15]. The beginning of the dashed arrows pointing up in (b) shows the beginning of the 
new phase. The existence of phonons in 7-phase in (b) are marked with the solid arrows pointing down. The a-phase extends up to 
820°C; (c) and (d) temperature dependence (RT-500°C) of phonon frequencies variation and FWHM for 72 (peak 1), 140 (peak 2), 
and 171(peak 3) cm^^, respectively. 



ruled out. However, the Raman frequencies before and 
after thermal cycling remain unchanged makes this un- 
likely. This fact favors films over bulk or powder samples 
for very high temperature studies in the future. Reach- 
ing the tetragonal and cubic phases extrapolated from 
the powder study of Arnold et al. 0] does not seem im- 
possible with thin films. 



B. Phonon anomalies near Tn: 



There are discrepancies in the literature regarding both 
the crystal structure of (001) BFO thin films e.g. with 
several reports claiming tetragonal , rhombohe- 

dral m, IIO], and monoclinic [3l|, HI] structure, and 
its phonons. Of particular interest regarding phonon- 
magnon coupling in BFO was the report [24 ] of a very 
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phase transition, we followed very closely the temper- 
ature dependence of few intense phonon modes i.e 72 
(Peak 1, we considered it as single peak for the sim- 
plicity), 140 (Peak 2), and 171 (Peak 3) citT^ from 
RT up to 500° C (Fig. [l]a,). In general, the change in 
phonon frequency band and width with temperature can 
be caused by several factors, such as anharmonic scat- 
tering, spin-phonon coupling, lattice expansion and/or 
contraction due to anharmonicity and/or magnetostric- 
tion effects, and phonon renormalization resulting from 
electron-phonon coupling [sj]. The latter one is not ap- 
plicable here as BFO is a highly resistive material, and 
the carrier concentration is low. The change in ionic 
binding energies with temperature also affects the change 
in phonon bandwidth in ionic compounds. However, this 
is not applicable here as BFO insulating. Figs. [TJ; and 
d reveal the fluctuation of phonon frequency and full- 
width at half maxima (FWHM) around 370°C (onset 
could be ~360°C), which happens to be the Tn of BFO. 
Near TNwe see small (1 or 2 cm~^) changes in both peak 
frequency and linewidth for three phonon modes at 72, 
140, and 171 cm~^. These satisfy the non-mean-field 
predictions of Nugroho et al [3^. This behavior could 
be the manifestation of phonon-magnon interaction that 
vanishes above Tn, 370° C. This weak interaction can be 
explained with the fact that the magnetic phase transi- 
tion is not accompanied by a structural phase change. 
The observation of a rather weak phonon-magnon inter- 
action is also consistent with the theoretical prediction of 
weak magnetization and ME coupling in BFO thin film 
by Ederer et al^. 



C. Phonon anomalies near 140-150K: 



FIG. 2: Room-temperature XRD patterns taken using CuKc 
(1.5406A) radiation (a) and Raman spectra (b) of (OOl)BFO 
film on SRO/STO before and after thermal cycle up to 1000° C. 
XRD pattern of STO is given for comparison. 



large (40 cm~^) change in the frequency of one long wave- 
length phonon branch near Tn. We emphasize in the 
present work that we see no such phenomenon (Fig. [T^) . 
Instead we see in FigsdJ; and d very small changes in 
frequency (1 or 2 cm~^) and linewidth of several polar 
modes, and we model them according to the non-mean- 
field theory of Nugroho et al [3^ . 

In general, there are three kinds of phonon anoma- 
lies observed near the phase transition temperatures: a 
sigmoidal S-shaped change in frequency (such as that re- 
ported by Haumont et al. [2J| but not found in our work); 
a step discontinuity; or a small "bump" (increase) that 
returns to the background level a few degrees above or 
below the transition temperature. 

In order to study the evolution of Raman signature 
around the antiferromagnetic-paramagnetic (AFM-PM) 



In order to study the cryogenic behavior of BFO thin 
films, we carried out scattering measurements (Fig. [3^) 
down to 81 K; no significant change in Raman spectra 
has been observed, indicating the RT structure remained 
unchanged down to 81K. However, a close observation 
shows phonon anomalies around 150K (Figs. [5)3 and c). 
This agrees with the observation of change in magnetic 
order at 150K by Pradhan et al. [s^, but the nature of 
this phase transition remains moot. 

Figs, and c show small anomalies in the frequency 
and linewidth of two phonon branches at 140 and 171 
cm~^near 140K, a temperature at which anomalies have 
previously been reported. Although the changes are 
small and the data sparse, they are highly reproducible. 
We are aware of similar observations on single crystals 
by Brahim Dkhil (37| , and we thank him for preprints of 
his work. The present data merely show that the same 
effects are present in thin films. 

The nature of the phase transition at 140K remains un- 
known. There are anomalies in magnon scattering cross- 
sections j3S ,.3iS and linewidth [40, 41], in mechanical loss 
tangent [3a, [i^l and 140K is the end-point in Almeida- 
Thouless data plots [40,|4l|. However, the earlier sugges- 
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tion by our group that 140K is a spin-reorientation 
transition temperature is not confirmed by very recent 
neutron scattering studies [13] and spin-glass effects [45j| 
have also been suggested but are unproven. 

D. Weak phonon-magnon coupling: 

In our Raman spectra (Figs. [T]: and d) we observed 
small (1 or 2 cm~^) increases or decreases in phonon 
frequencies very near Tn. Note that this is observed 
for several different phonon symmetries. The symmetry- 
independence of the phonon-magnon coupling implies an 
interaction of form P^M^ (where P and M are the po- 
larization and magnetization, respectively) in the free en- 
ergy, as first su gges ted for magnetoelectrics by Smolen- 
skii and Chupis . In general the coupling of phonons 
and magnons can occur through several different micro- 
scopic physical models: The Torrance-Slonczewski model 
[46| involves modulation of the crystal field at the spin 
site by the optical phonon and is significant for ions with 
unqucnched orbital angular momenta, such as Co^^ or 
Fe"*"^; the model of Buyers et al. is an angular mo- 
mentum coupling of spins in octahedra where the opti- 
cal phonon eigenvector is rotation-like as viewed from 
the magnetic ion. However, such models do not give 
frequency anomalies near Tn hke those we observed in 
BiFeOa . A rather detailed model of magneto-capacitance 
was given by Fox et al. [4811 for BaMnF4 near Tn, and 
related models by Scott [43 and by Glass et al. [50| for 
the BaMF4 family near T(2D), the two-dimensional spin 
ordering temperature [typically ca. 3Tn in that family]. 
The free energy of Fox et al. can be defined as: 

G = fiL^,Ll) + BM^ + (bo + hp + b2P^)M,L, + (1) 

where P is along y, the polar axis and z is the sublat- 
tice magnetization direction; L is {l/2)gfi S^/^ (S+1)^/^ 
[ XI (^j " down)] and M is the weak magnetization 
M = g^ (S+l)i/2 (Sj up + Sj down)]. Note 
that p is not the total polarization P but only the part 
induced by the magnetoelectric coupling: P = P,- + p 
and it was carried out to second order in polarization 
P, sublattice magnetization M, and weak ferromagneti- 
zation L, with the result that the magneto-capacitance 
varies with temperature as (bo - bib2)L2(T), where bp, 
bi, and b2 are respectively the coefficients of magneto- 
electric free energy terms independent of, linear in, and 
quadratic in polarization P. The authors noted that the 
mean-field theory, although generally not satisfactory for 
magnetic transitions, works well for weakly canted fer- 
romagnets because the expansion parameter L is small 
at all temperatures. Note that the sign of the magneto- 
capacitance term can be positive or negative depending 
upon the magnitude of (bo - bib2). Because they used 
mean field theory, their work neglected the small term 
near Tjq due to fluctuations considered below. 

Although the second-order theory of Fox et al. [1^ was 
satisfactory for describing all the data in BaMnF4 near 



Tn, it is not sufficient for phonon behavior in BiFeOa. In 
this case it is necessary to go to fourth order in L. The 
reasons are explained by Nugroho et al. [33| in their work 
on YbMnOa . In this case the key term in the free energy 
is of the form gP^L^, which for weak coupling gives an 
explicit interaction of electric field to L that results in a 
magnetocapacitance of 

(g^PykT) J [< L^ix)L^{0) >-<L^ >^]dx. (2) 

Although this fourth-order term is higher order than the 
terms in g<L^> considered by Fox et al. [il], it is sin- 
gular at Tn, because it is proportional to the cube of the 
correlation length (77) that diverges at Tn. 

The result is that the phonons in Raman effect in 
BiFeOa of any symmetry will be expected to have small 
anomalies in their frequencies at Tn. These small dips 
or jumps will be proportional to temperature t^°'~^^ [33] 
\5l\ . where t is reduced temperature, t = [TN-T]/TNand a 
is the critical exponent describing divergence of the spe- 
cific heat [5^ . Since the a is typically small, the phonon 
frequencies should vary approximately as Tn / [Tn-T] near 
Tn and in principle could be used to evaluate critical ex- 
ponent a. However in the present work the data are too 
imprecise for this chore, and even higher resolution would 
by insufficient due to phonon damping. The bump in 
phonon frequencies and linewidth are qualitatively pre- 
dicted from the non-mean field theory of Nugroho et al. 
[SSj . However, their model does not predict magnitudes 
for the height (increase in frequency) or width (how near 
the transition the increase occurs) of the bump. No 
anomaly at all is predicted by the mean-field theory of 
Fox et al. [i^, which does not consider terms in the 
free energy introduced by Nugroho et al. A similar 

behavior has been observed at 150 K (Figs. [5)3 and c) 
could be due to the change in magnetic ordering ^SQj . As 
matter of coincidence the bottom electrode SRO has a 
ferromagnetic phase transition at 150 K. Note that none 
of these peaks is related to SRO and a modulated effect 
is highly unlikely, but not impossible. 

In summary, our Raman frequencies near Tn and 150 K 
show small peaks or dips for all phonon modes that are 
qualitatively similar to those predicted by Nugroho et al. 
[33i] . implying a general interaction of form P^L^, and the 
need for a fluctuation term neglected in the mean-field, 
weak ferromagnetism model of Fox et al [i^ . 

IV. CONCLUSION 

In conclusion, high quality epitaxial (OOl)BFO films 
have been grown on (100) STO substrates using PLD. 
The XRD studies showed that films are c-axis oriented 
with high degree of crystallinity. The RT polarized 
Raman scattering of (OOl)BFO films showed pseudo- 
orthorhombic monoclinic crystal structure contrary to 
the rhombohedral and tetragonal symmetries reported 
earlier. We observed the ferroelectric phase transition 
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FIG. 3: (a) Temperature dependent Raman spectra of (OOl)BFO film on SRO/STO from 81K-RT; (b) and (c) temperature 
dependence of phonon frequencies variation and FWHIVl for 72 (peak 1), 140 (peak 2), and 171 (peak 3)cm^^ down to 81K from 
RT, respectively. 



at around 820° C and no softening of Raman modes was 
observed at low frequencies, as in BFO single crystals. 
The AFM-PM phase transition at around 370°C caused 
some small changes in the phonon frequencies, linewidth, 
and/or intensities of several low frequency modes, in- 
dicating ME coupling in the material. A non-cubic 7- 
BiFeOs phase was observed between 931-950°C in the 
BFO thin films, in agreement with the accepted BiFeOa 
phase diagram. The spectra also revealed an evidence of 
a cubic (5-phase around 1000° C in thin films that is not 
stable in powder and bulk. 
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